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Abstract 

We analyze the phenomena of condensate collapse, as described by Donley et al 0121, by focusing 
on the behavior of excitations or fluctuations above the condensate, as driven by the dynamics of 
the condensate, rather than the dynamics of the condensate alone or the kinetics of the atoms. The 
dynamics of the condensate squeezes and amplifies the quantum excitations, mixing the positive 
and negative frequency components of their wave functions thereby creating particles which appear 
as bursts and jets. By analyzing the changing amplitude and particle content of these excitations, 
our simple physical picture explains well the overall features of the collapse phenomena and provide 
excellent quantitative fits with experimental data on several aspects, such as the scaling behavior of 
the collapse time and the amount of particles in the jet. The predictions of the bursts at this level 
of approximation is less than satisfactory but may be improved on by including the backreaction 
of the excitations on the condensate. The mechanism behind the dominant effect ~ parametric 
amplification of vacuum fluctuations and freezing of modes outside of horizon - is similar to that of 
cosmological particle creation and structure formation in a rapid quench (which is fundamentally 
different from Hawking radiation in black holes). This shows that BEC dynamics is a promising 
venue for doing 'laboratory cosmology'. 
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1 Introduction 



We introduce a new perspective in the analysis of the phenomena of condensate collapse, described 
by Donley et al ^[2], by focusing on the behavior of fluctuations above the condensate, rather than 
the condensate itself. We show that the condensate dynamics squeezes, amplifies, and mixes positive 
and negative frequency components of the wave functions of the condensate excitations. In addition 
to providing a good qualitative understanding of the general picture our theory also produces precise 
predictions, specifically, on the critical number of particles at the first instance when the instability 
sets in, the scaling of the waiting time tcoUapse and the number of particles in a jet. In this rendition we 
point out the analogy between the evolution of quantum excitations of the collapsing condensate and 
the vacuum fluctuations parametrically amplified by the background spacetime in the Early Universe, 
suggesting a new venue for "laboratory cosmology" . 

A condensate formed from a gas of cold (3nK) Rubidium atoms is rendered unstable by a sudden 
inversion of the sign of the interaction between atoms. After a waiting time t^oUapse^ the condensate 
implodes, and a fraction of the condensate atoms are seen to oscillate within the magnetic trap which 
contains the gas (see below and ^). These atoms are said to belong to a "burst". In the experiments 
described by Donley et al. - to single out this controlled BEC collapse experiment from the others, 
we shall adhere to the namesake Bose Novae - the interaction is again suddenly turned off after a 
time Tevoive- For a certain range of values of Tevoive, new emissions of atoms from the condensate are 
observed, the so-called "jets". Jets are distinct from bursts: they are colder, weaker, and have a 
distinctive disk-like shape. 

The Donley et al. experiment takes full advantage of the tunability of the effective atomic interac- 
tion due to a Feshbach resonance characteristic of ^^Rubidium [3111 ■ The resonance is caused by the 
presence of a bound state whose binding energy may be tuned by means of an external magnetic field. 
In later experiments [11101) observed fluctuations in the number of particles in the condensate have 
been well-explained as arising from oscillations between the usual atomic condensate and a molecular 
state |71IHl[i[inilIIllI21. 

These oscillations were observed for magnetic fields in the order of 160G, where the effective 
scattering length is of the order of 500ao (and positive) (oq = 0.529 lO^^'^m is the Bohr radius) and 
the frequency of oscillations is of hundreds of KHz [SI [El- By contrast, in the Donley et al. experiment 

typical fields were around 167G, the scattering length was only tens of Bohr radii (and negative) 
and the frequency of atom - molecule oscillations may be estimated as well over ten MHz . Under 
these conditions it is unlikely that the molecular condensate plays any important dynamical role, and 
indeed no oscillations are reported in the original paper (for the opposite view, see Jll). For these 
reasons and to highlight the mechanism particular to this experiment, we shall not include explicitly 
a molecular condensate in our model but discuss in detail the one - field model. However, if need to, 
this may be done in a very simple way, by including a second field to describe molecular destruction 
and creation operators [HI IIUI ITlj . We will elaborate on this point in a later subsection 

There is a vast literature attempting to provide theoretical explanations of collapsing condensates 
^1 El El lEl El Uni [^ • In addition to speculations that Bose Novae is due to molecular oscillations 
as alluded above (which we view as of secondary importance) the most serious theoretical attempt 
is based on the Gross-Pitaevskii equation with explicitly introduced nonlinear terms to account for 
multiparticle interactions |2HI22j. We will show that the primary mechanism responsible for the main 
features of the Bose Novae experiment originates from the dynamics of quantum fluctuations around 
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the background condensate field(s). We start with the Heisenberg operator for the many body wave 
function and spht it into a c-number part describing the condensate amphtude and a q-number part 
describing cohective excitations (not individual atoms) above the condensate. We then derive an 
evolution equation for the wave function operator of the quantum (non-condensate) excitations under 
an improved Hartree approach, the so-called Popov approximation [2'6l I24L f7!jl \2ij\ . 

In this paper, we use a "test field" approximation, by adopting (rather than deriving) the specific 
evolution of the condensate extracted from the experiment as given and study the dynamics of the 
excitations riding on this dynamics. Note that the experimentally given condensate dynamics is 
different from the mean field dynamics obtained from a solution of the GP equation, because the 
former includes the dynamical effects of the fluctuations. Finding a self-consistent solution of the 
evolution equations for both the condensate and its fluctuations is called the 'backreaction problem'. 
It has been studied in detail in problems of similar contexts such as cosmological particle production 
(see below). Theoretical investigations for BEG fluctuations dynamics can be found in Refs. [71 1271 05] . 
The squeezing of quantum unstable modes and its back reactions on the condensate has been considered 
before, e.g., as a damping mechanism for coherent condensate oscillations j29J , and also applied to the 
description of condensate collapse j^,^ J2 , 33 , Field theory methods have recently been applied 
to the problem of formation and stability of Bose condensates |351 136| . Fluctuations have also been 
considered by Goral et. al. |37j and Graham et al. [HE]- Our work differs from them in the emphasis 
we place on the behavior of the quantum excitations as a consequence of condensate dynamics. 

Particularly relevant to the present work is Ref. where condensate collapse is analyzed from 
the point of view that the physics is mainly due to the dynamics of quantum fluctuations, the same 
view as we hold here. There, the trapping potential is replaced by a normalizing box, whose volume 
is eventually taken to inflnity. Our analysis in Sections 2 and 3 is for a more realistic geometry, 
which enables us to compare quantitatively to experiments. The analysis of bursts versus jets given 
in Section 4 however originates from a new concept inspired by cosmological processes. 

To the extent that many phenomena observed in connection to the collapse of this nature (Bose 
Novae) are essentially the result of a quantum fluctuation fleld (the non-condensate) interacting with 
a time dependent background (the condensate), as we believe it is, there is a close analog with similar 
processes in the early universe, speciflcally, vacuum particle creation from an time dependent external 
fleld ISni or in a curved background spacetime [lUlEI- (For a squeezed state depiction of this process, 
see, e.g., ^2j and references therein.) One could view condensate collapse as a laboratory realization 
of cosmological particle creation during quenching. (Note this is not the physical process behind 
black hole particle creation, as in the Hawking effect, much attention drawn to its detection in BEG 
notwithstanding US]-) In this process there is a competition between two (inverse) time scales, the 
physical frequency of the mode under consideration, and the inverse collapse (expansion) rate of 
the condensate. In cosmology the inverse expansion rate is the Hubble constant for the background 
spacetime. While a mode whose physical frequency is higher than the Hubble constant, we refer to 
it as "inside the horizon", and its behavior is oscillatory. When the converse obtains, the mode is 
"outside the horizon". They are depicted as 'frozen' because they do not oscillate (see below), but are 
amplifled 01] . This ampliflcation is largely responsible for the observed primordial density contrast 
in the Universe |45j . 

In the Bose-Novae collapse problem, the role of the "Hubble" constant is played by the inverse 
growth (exponential) rate of the most unstable mode of the condensate, which is determined by the 
instantaneous number of particles in the condensate. Modes whose natural frequency is greater than 
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the corresponding scale are relatively impervious to the dynamical condensate, but when the converse 
obtains, consequences are drastic. When the exponential growth is the dominant factor, the mode 
is frozen; instead of oscillating, it is being amplified, a process which is analogous to the growth of 
fluctuations during spinodal decomposition j46j . 

In the same way that modes that left the horizon during inflation return during the radiation and 
matter dominated eras, giving rise to acoustic oscillations, as the unstable condensate sheds its atoms 
and approaches stability, the band of "frozen" modes narrows: we say that modes "thaw" as they turn 
from exponentially increasing to oscillatory behavior. The crux of the matter is that only oscillating 
modes are detected through destructive absorption imaging (see below). Whenever a mode thaws, it 
is perceived as if particles were being created. In the conditions of the experiment the initial number 
of actual particles above the condensate is negligible, and we may describe the phenomenon as particle 
creation from the vacuum. 

To summarize, the key idea in our understanding of the phenomena associated with a condensate 
collapse is that of a dynamical background field of the condensate squeezing and mixing the positive 
and negative frequency components of its quantum excitations, thereby creating particles from the 
vacuum. The viewpoint of this work may be easily incorporated in a first principles approach as taken 
in e.g. |471 1481 1151 l5Uj . The remarkable analogy between condensate collapse and quantum processes 
in the Early Universe and spinodal decomposition in phase transitions may stimulate new related 
experiments in BEC to be carried out to address these problems in cosmology and condensed matter 
physics [5T| . 

This paper is organized as follows. Section II we briefly review the phenomenology of condensate 
collapse and set up the basic mathematical model. In Section III we give a discussion of the onset of 
instability and of the scaling of the waiting time tcoiiapse with the scattering length. In Section IV, we 
turn to a discussion of bursts and jets, based on the distinction between frozen and thawed modes. By 
postulating a specific condensate evolution (extracted from the experiment) we obtain quantitative 
predictions for the number of particles in a jet as a function of the time revolve (when the scattering 
length is brought to zero). Our results are summarized in the final Section. A few technical details 
are left to the Appendixes. 

2 The model 

2.1 An overview of condensate collapse 

In the Bose Novae experiment a gas of ^^Rb atoms at a temperature of 3 nK is prepared in 
a state where they behave essentially as a free gas within an anisotropic harmonic trap (see Figs. 1 
and 2). The trap has a cylindrical geometry [let r = {p^ip,z) be the usual cylindrical coordinates, 
with ip the azimuthal angle]; the trap frequencies are Vradiai — 

17.5 Hz {ijJradial = '2'Tn^radial = HO 

Hz) and Vaxiai = 6.8 Hz {uaxiai = 42.7 Hz). At time t = 0, the scattering length a (see below) is 
suddenly turned to a negative value —acoiiapse- This configuration is known to be unstable whenever 
the number Nq of atoms in the condensate exceeds Kaho/\a\, where Uho is a characteristic length of 
the trap. The coefficient k was reported as k = 0.46 [HS], but later measurements suggest it should be 
raised to k = 0.55 

For Nq = 6000 atoms, the instability threshold is reported at |a| = Ucr = 5.12 oq, where oq = 0.529 
10^^'^m. Therefore aho = 5.6 IO^oq = 3 10~^m. If we write a^^ = h/Mco and introduce the atomic mass 
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of 85mp ~ 1.3 10 ^^kg and h = 1.05 10 ^^Js, this corresponds to a frequency u = 90Hz. This agrees 
well with the geometric average u = {^axial^radial) ~ 80Hz. It follows that t^o = uj~^ = 12.5ms. 

In spite of the instability, over a time tcoiiapse there is no significant decay of the condensate. 
tcoiiapse depends very strongly on acoiiapse (see below). After tcoiiapse^ the number of atoms in the 
condensate falls exponentially. If left to itself, the condensate eventually stabilizes retaining a number 

Nremnant of atomS. 

During this period, a cloud or burst of atoms is observed. These atoms oscillate in the radial 
direction with the trap radial frequency. The energy associated with the burst is larger in the radial 
direction than in the axial direction. The number of atoms in the burst increases with the time elapsed 
since the decay begins, with a time constant of about 1.2 ms; after 7 ms, the number of burst atoms 
stabilizes. About a fifth of the atoms in the condensate go into the burst, with variations of about 
20%. 

Condensate decay is interrupted at a time t = T^volve^ when the scattering length is again tuned 
to a positive value. If the condensate is already stabilized at Tf,yoive-, nothing too drastic happens, but 
otherwise a new phenomenon appears, namely the expulsion of a jet of atoms. Jets have much lower 
kinetic energy than bursts (a few nano-Kelvin) . 

2.2 Basic equations^ 

The model is based on the Hamiltonian operator for N interacting atoms with mass M in a trap 

with the total number operator N given by 

N = j dv (2) 

Here V (r) is the trap potential and U is (assumed to be a short ranged) the interaction between the 
atoms. We introduce a dimensionless field operator ^ (x) 

^{r)^a-i^'^^{x) (3) 
where aho is a characteristic length of the trap, and a dimensionless coupling constant u 

U = hcoal^u (4) 

In terms of the scattering length a (which we define as positive for an attractive interaction), we 
have 

(^) 

""^In writing down tlic main equations of our model, we choose a sign convention which makes the effective couphng 
constant positive for an attractive interaction, and a system of units adapted to the problem. Taking the average frequency 
u} as reference, we may define length and time scales auo and tho (see above) and an energy scale E^o = Jiuj = MiJ^a\g. 
Prom now on, we shall choose units such that these three scales take the value 1. 
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So 



u = 



47r 



aho 



(6) 



The Hamiltonian and the trap potential can also be written in terms of dimensionless variables 

U=EhoH, Y{v) = EhoV{r) (7) 
Assuming a cylindrical shaped potential 

1 



V {v) = -{uy + u^y) 



(8) 



with radial p and longitudinal z coordinates measured in units of aho, with associated (dimensionless) 
frequencies uj^ = ^axiall'^ ~ 1/2 and = u^radiall'^ ~ a/2- 
* obeys the equation of motion 



and the equal time commutation relations (ETCR) 



whereby 



Here 



dt 



H = --VyV{r) 



(9) 



(10) 



(11) 



(12) 



is the (dimensionless) one-particle trap Hamiltonian. We decompose the Heisenberg operator ^ into 
a c-number condensate amplitude $ and a q-number noncondensate amplitude ■0, consisting of the 
fluctuations or excitations 



Then 



and 



*(r,i) = $(r,i) +V'(r,i) 
= |$|2 ^ ^^^t + ^) + 



(13) 
(14) 



n =<*%>= |$|^+ < >= "0 + (15) 

where () denotes expectation value with respect to the initial quantum state. The cross terms vanish 
because the mean field is defined with zero fluctuations. Here no = |$|^ is the number density of the 
condensate and h = is the number density corresponding to the fluctuations. Further, 
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2,At 



|$|2 ^ + 2 |$|^ V + + 2$VV + + V'V 



Under the self-consistent mean field approximation |^ 



n 



where fh is sometimes called the anomalous density, and 



(16) 



(17) 



(18) 



(19) 



Taking the expectation value of the Heisenberg equation we obtain the Gross - Pitaievsky equation 
(GPE) 



9$ 9 

i— - + u |$| $ + n{2$n + <^*m] = 



(20) 



Subtracting this from the original equation Eq. H1U() . and neglecting nonlinear terms we obtain the 
equation for the fluctuation field 



i— - H + 2um^ 
ot 



tp (x) + uf^'^tp'^ {x) = 



(21) 



This is the basic equation we shall use here for analyzing the behavior of the excitations (fluctuations) 
for a given condensate evolution. From the evolution of the fluctuations we can then calculate the 
modified evolution of the condensate. Note that only in taking the Bogoliubov approximation (setting 
n = m = 0) will one get a closed equation, the Gross-Pitaevskii equation, for the condensate. In 
principle these two equations for the fluctuations (PT|) and the condensate PU)) need be solved together 
in a self-consistent manner. (This is called the 'backreaction problem' first studied in external field 
or cosmological particle creation problems). Here we solve only the equation for the fluctuation field 
using the experimentally measured value for the background field. As mentioned earlier this is already 
a (lowest order) backreaction- modified or self-consistent solution which is an improvement over the 
mean field solution obtained from the GP equation. 
We next parametrize the wave functions as 



(22) 



ip = ipQe 



-ie 



(23) 



where and are real. We seek to express the two equations for $,■0 by the three equations for 
the three real quantities ^o,ipo,Q. Observe that 



(24) 
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therefore, 



^ + 1 (Ve)^ - -^H<^o + u^l + u {2n + Re m} = 
ot 2 <Po 

^ + VeV«>o + - (V^e) $0 + u^olm m = 
ot 2 V / 



The third equation may be simphfied by using the first 

We have also neglected terms which were nonlinear in fluctuations. 



iPo + u^I^pI = 



(25) 
(26) 
(27) 

(28) 



3 Onset of collapse and scaling of tcoiiapse 
3.1 The early stage 

During the first few miliseconds of evolution, we may regard the condensate density as time indepen- 
dent, and the condensate phase as homogeneous (see Appendix A for justification). We may then 
write the equation for the fluctuation field Eq. 1)21^ as 



d 

i--H + Eo 
ot 



1 







To solve this equation we decompose it into a self-adjoint and an anti-adjoint part 



with each part satisfying an equation 



dr] 
'di 



+ 



H-Eo-2u% C = 0. 



In the subspace orthogonal to the condensate, the operator has an inverse, whence 
and since the trap Hamiltonian is time - independent, we have 



dt^ 



(29) 
(30) 

(31) 
(32) 
(33) 

(34) 
(35) 
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Here 



Heff = H-Eo-2u^l 



(36) 



is the Hamiltonian for a particle moving in the potential 





n 



Then (t) will be a superposition of two harmonic oscillations with frequencies ±y/X^. To have an 
unstable condensate it is necessary that at least one of the is negative; the boundary of stability 
occurs when the lowest , say Ag is exactly zero. Since [H — Eq] has an inverse, the equation 
[H - Eq] Hefftpo = implies -ffe//V'o = 0. 

In conclusion, the onset of collapse corresponds to the minimum value of the condensate density 
for which the effective Hamiltonian Hgff develops a zero mode. 

3.2 The onset of collapse 

Wc shall now study the spectrum of the operator H^jf. The idea is that the low lying states, which are 
the most relevant ones to our discussion, will try to keep close to where the potential is a minimum, 
namely, the origin. If the width of the state is small enough, the potential will be nearly harmonic. 

One further consideration is that we are interested in the part of the fluctuation field which remains 
orthogonal to the condensate, since fluctuations along the condensate mode may be interpreted as 
condensate fluctuations rather than particle loss. This is granted for all modes with odd parity; for 
modes of even parity, it means that a certain combination must be excluded. The ground state of 
Hf,ff is certainly not orthogonal to the condensate, since neither have nodes. 

A simple but remarkably accurate way of determining whether -ffe// admits a zero mode is to 
consider its expectation value {Hf,ff) with respect to some appropriate trial state. The first excited 
eigenvalue -Eioo is the minimum possible expectation value, taken with respect to any normalized state 
orthogonal to the ground state. Therefore, if we may exhibit a wave function leading to a negative 
(Heff) , then necessarily E'loo is negative as well. 



Let us write 



H^ff = H§ + HP - 2u^l - Eo 



(39) 




(40) 



(41) 
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Recall for the condensate 



$2 ^ ^ -[....2+^pp2] 

^3/2'' 



(42) 



We shall try to minimize the expectation value of H with respect to a wave function of the form 



V^^^^ (z) Vo ^ (x) V'o" (y) 



(43) 



where iIjq' and "01^ are the fundamental and first excited states of an one-dimensional harmonic oscillator 
with arbitrary frequency O 



4J7\V4 



vr 



The free part yields 



4 I ' 0, 



The condensate part yields 



where 



(44) 
(45) 

(46) 
(47) 

(48) 

(49) 
(50) 

If we adopt the values lOz = 1/2, = \/2; relevant to the JILA experiment, then we obtain 

A a . Q 7g ^ Y.l. From the definition of A, 



d\ (r) [V'l^^ (z)# {x)^\,f (y) 
Putting them all together 



7r3/2 



(Heff) = 2 



■ to. 



A 



A = iu 



No 
.3/2 



TT' 



16 a \ 

0/10/ 



{Hf,ff) = for the first time when A = Acrit ~ 4.6, Q 
we conclude that instability will occur when 



O'crit 



A 



crit 



0.51 



(51) 

aho 16 

This result compares remarkably well with the experimental value k = 0.55, as well as with the 
theoretical estimate presented in Ref. j53| . This agreement may be seen as natural, as the equations 
we postulate for the fluctuations may be obtained from the linearization of the GPE, discarding both 
h and hi. In both calculations, the geometry of the trap plays a fundamental role. 
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3.3 Scaling of tcoiiapse 

As we have already noted, even for condensate densities above the stabihty hmit, no particles is seen 
to be lost from the condensate during a waiting time tcoiiapse- 

Even in the absence of a detailed model of the condensate evolution, the above analysis allows us 
to make a definite prediction of the way tcoiiapse scales with the scattering length. The basic idea is 
that, while tcoiiapse depends in a complex way on several time scales, some intrinsic to the condensate 
and some related to the condensate - noncondensate interaction, up to tcoiiapse the time scales intrinsic 
to the condensate are very large compared to the other processes. The nontrivial point is that even 
in this limit the time scales of the noncondensate remain finite, and so they fix the scale for tcoiiapse 
itself. Using the exponential growth of the first excited state as a measure, we are led to the estimate 

tcoiiapse ~ ^ 

Consider a value of A close to, yet higher than, the critical one. The fastest growing mode is 
the eigenvector of [H — Eq] Heff with most negative eigenvalue, say —a. When A takes the critical 
value Acr, cr = 0. Close to the critical point, we may compute a using time independent perturbation 
theory, even if [H — Eq] Hcff is not Hermitian. The idea is that at the critical value of A there is a 
single normalized state |100) orthogonal to the condensate such that [H — Eq] Hcff |100) = and 

(100| [H - Eq] Hcff m + (V'l [H - Eq] Hcff |100) = (52) 

for any state orthogonal to both the condensate and |100) (otherwise, there would be states 
orthogonal to the condensate with {[H — EQ]Hcff) < 0, which is impossible). Now for A slightly 
above the critical value, let \a) = |100) + \6a) , (100| 6a) = be the eigenstate corresponding to the 
minimum eigenvalue —a. Then 

-a = (100| [H - Eq] Hcff \a) = (100| [H - Eq] Hcff \100) + higher order (53) 

and so 

(^-l)([^-^o]2ufg) (54) 

where the expectation value is computed at the critical point. Since we know the form of the wave 
function, computing the expectation value is a simple exercise. We find that the growing mode 
increases as expet, where 

/A 

and T = ([H — Eq] 2u^q) ~ 0.6. Since A/Acrit = o/ocr for the given total number of atoms, 

tcoiiapse ~ tcrit I 1 ) (^6) 



cr 



The power law Eq. (|56|) describes with great accuracy the way tcoiiapse scales with the scattering 
length; the actual prediction tcrit ~ is correct only as an order of magnitude estimate. In natural 
units, is about 20ms, while the best fit to the experimental data is obtained for tcrit ~ 5ms. 
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We should stress that the factor of 4 discrepancy between our prediction and the experimental 
result could be easily explained away in terms of a more complete model of condensate depletion. For 
example, the probability for two particles in the condensate to collide and transfer to the growing 

mode, Bose enhanced by the population in the latter, would go as ^V'^V') ) thereby yielding the factor 
of four in the exponent. However, in this paper we shall keep to our goal of providing the simplest 
theory which gives the best qualitative description of the Bose Novae phenomena, hereby understood 
as arising from the quantum dynamics of the fluctuations. 

An analysis of collapse based on the Popov approximation equations, for a simplified geometry, is 
given in However, the time i^vL ~ (^^-^o)~^ .54 , which is herewith identified as characteristic of 
collapse, does not account for the enhancement of t collapse near the critical point. 

In Fig. 3 we plot the scaling law H56|) (full line) derived here and compare it with the experimental 
data for A^o = 6000 as reported in Refs. 012] (small black points), the tNL ~ iuNo)~^ prediction 
(suitable scaled) as given in [2] (dashed line) and the results of numerical simulations reported 
in |2j (large grey dots). While all three theoretical predictions may be considered satisfactory, the 
tNL ~ (wAq)"^ fails to describe the divergence of tcoUapse as the critical point is approached. The 
results of numerical simulations reported in |2j| based on an improved Gross-Pitaevskii equation tend 
to be systematically above the experimental results. In a classical instability, the unstable modes 
must grow from zero, while in a quantum instability, they are always seeded by their own zero - point 
fluctuations, which speeds up the development of the instability. Therefore, the fact that numerical 
simulations based on classical instability tend to overestimate tcoUapse may be a further indication of 
the quantum origin of the phenomenon. 

3.4 Minor role of a molecular condensate 

Having obtained some insight into the early stages of collapse, let us discuss briefly how the present 
model may be modified to explicitly account for the population of a molecular state. We refer the 
reader to Refs. ^[EHEH for further details. 

The system of atoms and molecules at the Feshbach resonant state may be described by a theory 
where the fundamental degrees of freedom are an atomic field ^ (x) (we use the dimensionless am- 
plitude introduced in the previous Section) and a molecular field A (x) . The Hamiltonian takes the 
form 

H = Hatom + Hmol + -f^fee + Hint (57) 

where 

Hatom = Htrap + H^kgd- (58) 

Here Htrap given in Eq. ((T^ describes the dynamics of free atoms in the trap, Hbkgd accounts for the 
background interaction between atoms (that is, very far from resonance) 



Hbkgd 



-Ubkgd\ 



j j (fx ^+2^2 (59) 



2 

where j6j (recall that our sign convention is a positive scattering length for an attractive interaction) 
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Ubkgd = 47r^^, abkgd = 450 Oq =^ Ubkgd = 0.1 (60) 

Hraol describes the motion of molecules in the trap. We treat a molecule as a particle of mass 2M 
(where M is the mass of the atom) subject to a potential 2V (where V is the trap potential seen by 
the atoms) [HS]! so the frequency of oscillations in the trap is the same for atoms and molecules. H^e 
accounts for the energy difference between a molecule and two atoms 

Hbe = -£ {B) J d^x A^A (x) (61) 

Plots of the binding energy e (B) as a function of the applied field B are given in |H1E]. Finally, Hint 
accounts for the formation and dissociation of molecules 



H.. 



int 



'J (a^^"^ + ^^"^a) (x) (62) 



The binding energy e (B) is the crucial input parameter in the model, e (B) = at the resonance 
field B = Bpeak = 155G jl3j and increases with larger S's. At values of i? ~ 156G typical of 
the Ramsey fringes experiment jH] we already have binding energies of the order of lOKHz; this is 
much larger than the lOOHz typical of motion in the trap, and therefore H^oi is negligible (in the 
condensate collapse experiment P typical fields where B ~ 167G, and e {B) was much larger). Under 
this approximation, the full Heisenberg equation of motion for the molecular field 

i^^A = -e{B)A + T^'^ (63) 

may be solved analytically 

A (r, t) = A (r, 0) e^^* - iV /* dti e*^(*-*')^2 q (g4) 

Jo 

// the atomic field is slowly varying in the scale of then the integral is dominated by the upper 
limit, and simplifies to 

A{r,t) = Afree{r,t) + ^^^{r,t) (65) 

The free part decouples from the atomic condensate, and the effect of the driven part is to introduce 
an effective interaction among atoms, so that u^kgd in H^^gd is replaced by 

Ueff = Ubkgd - (66) 

which is the familiar pattern which makes it possible to tune the scattering length. The effective 
scattering length vanishes at B = B^ero = 165. 75G ^Hl (see Fig. 1). 

We conclude that under this approximation, the number of molecules - disregarding those already 
present independently of the atomic condensate - is 
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If Natom ~ 16000 and Ueff ~ —Ubkgd inii then the molecular condensate becomes important when 
£ (B) ~ Ubkgd^atom^ho = 128kHz. This regime was achieved in p|. Under these conditions we cannot 
rely on the above approximation, but rather we must solve the model. This yields atom - molecule 
oscillations with a frequency set by £ {B) . 

In contrast, in the condensate collapse experiment the magnetic field was driven above B^ero: 
with typical values B = 167G. For these fields, the molecular binding energy is close to 12MHz 
(see jlHl EE]), or 15 10^ ujq Therefore Nmoi ~ 0-005Natom ~ 80. This is less than the number of 
atoms even in the earliest measured jets (see [I] and below), suggesting that no major role is played 
by molecular recombination. This conclusion is consistent with the fact that no oscillations in the 
number of particles in the condensate are seen (see and below). It is interesting to observe that the 
numerical simulations presented in Ref. ^] assume a set of parameters chosen to display the effect 
most clearly, not attempting to be realistic. 

This leaves open the possibility that the approximation may break down because the atomic field 
is not slowly varying with respect to £~^. For example, the time it takes to change the magnetic field 
from Bzero to Bf.^oive is probably of the order of 10^s= 8 10~^u;^^^ This is a short time compared 
to under the conditions of the Ramsey fringes experiment [HI, but a rather easy stroll in the 
condensate collapse experiment ^ . The next scale in which we expect the condensate to react (which 
is also of the order of the inverse atomic chemical potential) is t a^l |54| I31j , but this is of the order of 
milliseconds and therefore way too long. 

The conclusion is that formation or dissociation of molecules, under the conditions of the con- 
densate collapse experiment is unlikely to have played a central role. Moreover, in models where 
the molecular state is explicitly included, it is still necessary to account for the quantum dynamics 
of the atomic field. This is done, for example, in Ref. ^Ij by tracking the fluctuation two-point 
functions along with the condensates, or in by providing independent equations for the anomalous 
atom pair amplitude. Resonance of molecular field with the condensate alone, though remarkably 
successful in the later experiments, cannot account for the phenomena of this experiment. For the 
Bose Novae experiment, the dominant qualitative features are determined by the dynamical vacuum 
squeezing mechanism proposed here; including a molecular field can only provide a minor quantitative 
improvement. 

4 Evolution of fluctuations: Bursts and Jets 

In the previous Section, we were mostly concerned with identifying the factors which can render the 
condensate unstable. In this Section, we shall consider the quantum evolution of fluctuations, as a test 
field riding on the collapsing condensate extracted from experiment. The basic dynamical equations 
are, from before, Eqs. (jS3), p2|) and (|33|) . where now <I>q is assumed to have both space and time 
dependence. The initial state is defined by the condition that u = for t < 0; we shall take it to 
be the particle vacuum |0) , which is defined by ip{x,{)) |0) = everywhere. Rather than seeking a 
rigorous solution, we wish to understand the structure of the problem in order to display the extent to 
which the phenomenology of condensate collapse is determined by fundamental (yet simple) quantum 
dynamical effects such as parametric amplification and particle creation. With this goal in mind, we 
shall make some physically motivated simplifications. 

We have already seen that it is possible to introduce a mode decomposition of the ^ operator based 
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on the eigenfunctions of [H — Eq] H^f / . This shows that the preconceived notion that fluctuations wih 
react only to the local state of the condensate is flawed; the relevant modes are nonlocal and can 
sample conditions over a large portion of the trap. 

For short wavelengths A, we expect these eigenfunctions will approach trap eigenmodes, since 
H ~ A^^ >> 2u^Q. Also the fact that particles in bursts are seen to oscillate with the trap frequencies 
P suggests that their dynamics is determined by the trap Hamiltonian (see below). Of course, the 
trap eigenmodes would also be eigenfunctions of -ffe// if this operator and H commuted. Now, since 



we see that we may disregard the conmutator when the condensate is slowly varying with respect to 
the relevant modes. Early enough in the collapse, the typical scale for the condensate is aho, and 
this condition will hold for almost every mode. Thus we shall make the approximation of assuming a 
homogeneous condensate. The amplitude of the homogeneous condensate will be chosen to enforce the 
onset of instability. Since we expect the lowest trap mode to become unstable when = N^r = n/a 
(in units where aho = 1), we approximate 



where Nq (t) is the instantaneous total number of particles in the condensate. For constant Nq, this 
yields e ~ (a/acr — 1)^''^ for the time constant of the growing mode, which is close to the more 
rigorous estimate Eq. (|^^ above. 

In practice, in Ea. ()69() is a measure of the overlap between the condensate and the excitation 
modes, as in Eq. (|48|) above. Therefore, the approximation may be improved by adjusting k according 
to the range of modes where it will be applied. We shall adopt this practice in what follows. 

4.1 Evolution of fluctuations up to tcoiiapse 

We now proceed with the quantitative analysis of Bose Novae. We assume a given evolution extracted 
from experiment for the number of particles in the condensate and analyze the evolution of fluctuations 
treated test field on this dynamic background. 

Concretely, we shall consider evolutions where the number of particles in the condensate remains 
constant from the time t = when the scattering length is switched to a negative value up to some 
time tcoUapsei and decays exponentially from then on (see Fig. 4). In the actual experiments the 
condensate does not evaporate completely, but the number of particles in the remnant is much smaller 
than the initial number of atoms in the condensate, and is therefore negligible at the early times when 
the approximation of a homogeneous condensate is valid. 

Let us begin by considering the evolution up to the time of collapse tcoiiapse-, or the waiting period. 
Let Nq be the initial number of particles in the condensate, and Ocr = i^/Nq the corresponding critical 
scattering length. Trap eigenfunctions are labelled by a string of quantum numbers n = {nz,nx,ny) . 
The eigenvalues of the trap Hamiltonian are (with the zero energy already subtracted) Ef^ = nuj where 
id = {lOz,^^p, ^p) and nCo = tOz^z + ojp {rix + riy). We choose the eigenfunctions V'n to be real. 

Let us expand 




(68) 




(69) 
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Then 



^ = Z] «ff (t) 'il^n (r) 
n 

^ = Z] C^J (t) V'n 

n 

V = '^bft (t) ipfi 

n 

Cn = ^ («n + 4) 

— f Oi? — at) 
2i\ "/ 



d 
dt 







En - 








Enbn 



(70) 
(71) 
(72) 

(73) 
(74) 
(75) 
(76) 



We see that there are two kinds of modes, stable (oscillatory, or thawed) modes if En > (^) "^z) and 
unstable (growing, or frozen) modes if not. We shall consider each kind in turn. 

4.2 Stable modes and bursts 

First consider the case when En > (ir^) coz- We get 

Cfi{t) = ^ (an + at) (0) cos unt + (an - aVj (0) sinunt 



bn {t) = ^ ~ 4) (0) cos Lvnt - ^Xn (an + 4) (0) sinu;^jt 



where 



UJfi 









h 


En — 


( — ]oj^ 






\acrJ . 



XfiEn 



Xn = \h-- 



a \ iOz 
acr J En 



So 



where 



an{t) = fnit)an{0)+gnit)aUo), 



(77) 
(78) 

(79) 
(80) 

(81) 
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fn{t) = cos Lont + 



-C/r? 



sin ujfit 
2iXn 



9n{t) = -{ — 



Observe that 



(82) 
(83) 



l/n(i)l'-bn(i)l' = cos2a;^jt + 



Ctcr 



a iOz 

O'er 



sin ijjfit 



(84) 



as required by the commutation relations. Although we assume vacuum initial conditions, these modes 
do not remain empty. The density 



n{r,t) = E€M(4W«nW 



sin ojf^t 



(85) 



We see that the density has a constant term and an oscillatory term. In our view this oscillatory 
term is responsible for the appearance of "bursts" of particles oscillating within the trap. This point 
is worthy of some elaboration, as it will make clearer the contrast with the unstable modes to be 
discussed below. 

Let us approximate the amplitudes (r) by their WKB forms 



1=1-3 



V'i in) 



cos 



TT 

Di — Tlj — 



(86) 
(87) 



where 



Si 



dx Pi {x) 



Pi 



2niUJi - ajhf 



(89) 



After expanding all trigonometric functions, we see that the oscillatory part is a linear combination 
of terms of the form 



exp|2i [E 'S'i - ^^n^] } 



(90) 



in all sign combinations. If we look at the density at a given point and time (r, t) , we see that the 
modes which contribute effectively are those for which (in the limit Xn ^ 1) -^n S> ^z) 



iOi 



' ri dx 

Jo Pi (x) 



-t 



= 



(91) 
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This equation describes a particle which moves along a classical trajectory in the trap potential, 
starting from the origin at t = with momentum pi (0). If we include the phase shift in computing 
the saddle point, there is a time delay, but it is the same for all particles. 

We conclude that the oscillatory part of the density describes a swarm of particles moving along 
classical trajectories in the trap potential. These trajectories return to the origin at multiples of the 
inverse trap frequencies, thus producing a surge in the local density, which becomes large enough to 
be seen by destructive absorption imaging. These particles constitute the so-called bursts observed 
in the Bose - Nova experiment. 

4.3 Unstable modes and jets 

Let us now consider the opposite case Efi < ( ) ujz ■ Formally we may obtain the relevant formulae 
from the last section with the replacement Xft = i'&Hi = if^ni thus transforming cosoj^t coshu^f 
and sin w^t i sinh cr^t. We get 

Cn (t) = ^ {an + 4) (0) cosh cr^t + (a^ - at) (0) sinhu^t (92) 
bn (t) = ^ (on - at) (0) cosh afit + ^^9^ (o^ + at) (0) sinh a^i (93) 



So again 



But now 



On (t) = (t) an (0) + 5n (t) ai (0) (94) 



fn (t) = cosh afit + (l - I?!) sinh a^t (95) 

9H (t) = -^(l+^l) sinha^t = (—) ^ sinhcT^t (96) 

Physically the difference is huge. In the first place, the density is growing exponentially. But 
unlike the previous case, there is no oscillatory component. While the actual number of particles is 
increasing, there are no surges in the density caused by the sudden constructive interference of many 
modes. In particular these particles do not contribute to the central peak in the density distribution. 
In this sense, they can not be seen by destructive absorption imaging. Because these particles do not 
oscillate in the trap, in the sense above, we say these modes are frozen in the same sense used in 
theories of cosmological structure formation, i.e., that fluctuations in an evolving universe are said to 
freeze upon leaving the horizon |44j . 

However, these modes come alive at Tevoive-, when the scattering length is set to zero. Now they 
become ordinary trap modes, and oscillate in the trap in the same way as the burst described above. 
To the observer, they appear as a new injection of particles from the core of the condensate, which 
makes up the so-called 'jets'. The sudden activation of a frozen mode by turning off the particle - 
particle interaction may be described "thaw" . 

Observe that in this picture several conspicuous features of jets become obvious. Jets may only 
appear the turn - off time Tevoive is earlier than the formation of the remnant. Once the condensate 
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is stable again, there are no more frozen modes to thaw. On the other hand, jets will appear (as 
observed) for T^voive < tcoiiapse, when the condensate has not yet shed any particles. Also jets must be 
less energetic than bursts, since they are composed of lower modes. Because of their relative weakness, 
treating them as test particles riding on the dynamic condensate as is the approximation adopted here 
('test field') works better for jets than for bursts. A more accurate depiction of the dynamics of bursts 
may require consideration of the backreaction of the fluctuations on the condensate. 

4.4 Beyond tcoiiapse ■ thawing 

The main physics we presented in the last Section is that the amplitude for a normal mode in the 
excitation (fluctuations) of the condensate evolves as a harmonic oscillator with a complex frequency 
dictated by the dynamics of the condensate. Parametric amplification (or, in a quantum optics, 
squeezing) of the fluctuations populates the modes above the condensate, even from a vacuum state 
originally, manifesting as bursts or jets. Thus in our view the salient features of Bose - Novae, the 
specific phenomena arising from the collapse of a BEG formed from a dilute Bose gas arises from the 
squeezing of (vacuum) fluctuations by the condensate dynamics. 

Let us now discuss the behavior of fluctuations beyond tcoiiapse, when the number of particles in the 
condensate, and therefore the instantaneous frequency of the excited modes, become time dependent. 
As in the previous Section, we shall assume nevertheless that the condensate remains homogeneous, 
thereby confining ourselves to the early stages of collapse. 

For concreteness, let us assume that the number of particles in the condensate remains con- 
stant A'^o = iVo = 16, 000 up to t = tcoiiapse = 3ms, and then decays exponentially Nq (t) = 
N^exp (— (t — tcoiiapse) /t), with r = 6ms (see Fig. 4). The mean frequency of the trap ui = 80Hz 
corresponds to a temperature huj/ks = 52.5nK. The actual trap frequencies reported in Ref. P 
^radial = HO Hz and ojaxial = 42.7 Hz Correspond to temperatures 28. 012 (axial) and 72. 162 (radial) 
nK. These are relatively high with respect to the sample temperature of 3 nK. For this reason, the 
initial number of particles above the condensate is negligible, and we may assume that we are dealing 
with particle creation from effectively the vacuum. This is the rationale behind calling this process 
'squeezing of the vacuum'. We shall assume the scattering length is brought to a = 36ao. We shall 
approximate the effect of the condensate on the fiuctuations by a constant level shift as in Eq. (|69j) . 
The best fit to the experimental data is obtained for k, = 0.46, which is satisfactorily close to the 
experimental value of 0.55. 

Shifting the origin of time to tcoiiapse for simplicity, we write A^o {t) = -/Voexp (— t/r) . After expand- 
ing in trap eigenmodes as in the previous Section, we obtain the equations 




(97) 




d 



(98) 



Therefore 




(99) 
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This equation clearly displays the two kinds of behavior described above (an exact solution is provided 
in Appendix B). If Efi > (^^) , the mode is always oscillatory. If E^^ < (^^) , the mode is frozen 
at tcoiiapse-, but thaws when exp (— t/r) ~ Ef^a/auJz (see Fig. 5). During the frozen period, the modes 
are amplified, but they only contribute to bursts after thaw. If the evolution is interrupted while still 
frozen, they appear as a jet. 

We therefore conclude that the number of particles Nj^t in a jet at time T^^^oive is essentially the 
total number of particles in all frozen modes at that time. If we write as before 

«n (t) = fn (t) an (0) + QH (t) (0) (100) 

then 

Njet{t)= E 

This is plotted in Fig 6, from the exact solution in Appendix B, with the parameters given above, 
and compared to the corresponding results as reported in pQ. We see that the agreement is excellent 
at early times (up to about 6ms). For latter times, the model overestimates the jet population. This 
is due to the fact that, by not considering the shrinking of the condensate, we are overestimating 
the overlap between the condensate and the fluctuations, thus delaying the thaw. It nevertheless 
reproduces the overall slope of particle number with Tevoive-, which is quite remarkable considering the 
simplicity of the model. 

The see-saw pattern follows from the discreteness of the modes. Modes thaw at discrete times 
t^. Between one and the next, occupation numbers of the unstable modes are increasing, and so the 
number of particles in the jet is a growing function of Tevoive- When the next stabilization occurs, the 
particles in the now stable mode no longer contribute to later jets, and the particle number in the jet 
decreases by a finite amount. This pattern accounts for the fact that a jet from a later Tevoive may 
be stronger than earlier ones, and also for the large variation in the number of particles in jets with 
similar values of Tevoive- It should also be remembered that we are computing the expected number of 
particles, but that, in the highly squeezed state which results from the frozen period, the fluctuations 
in particle number are comparable to the mean number itself. 

5 Conclusion 

In this paper, we have applied insights from the quantum field theory of particle creation and structure 
formation in cosmological spacetimes and the theory of second order phase transitions to a specific 
scenario of controlled collapse of a Bose-Einstein condensate, the so-called Bose Novae phenomena. We 
have described these phenomena as resulting from particle creation from the vacuum, induced by the 
time dependent condensate. This time dependence squeezes and amplifies the field operator describing 
excitations above the condensate. A key concept in our analysis borrowed from theories of cosmological 
structure formation is the drastic difference in the physical effects of frozen versus oscillatory modes: 
those whose physical frequencies are higher than the collapse rate oscillate and are rather impervious 
to the condensate, while those below (frozen modes) grow in time and get amplified, in a way similar 
to the growth of fluctuations during spinodal decomposition. As the condensate stabilizes and the 
collapse rate decreases the frozen modes begin to thaw. The appearance of oscillatory modes (in 



9n (ty 



(101) 
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second quantized language) is described as particle creation appearing in jets and bursts, as described 

in detail above. 

In order to focus on the key ideas we have adopted a number of simplifying assumptions. We 
take the condensate evolution as a given input from the experiments, rather than deriving it from 
fully self-consistent equations. We have treated excitations within the Popov approximation, which 
improves on the Hartree approach but is known to break down as the number of particles above the 
condensate increases. We have neglected the coupling between different excitation modes, considering 
only the coupling of each to the condensate. 

These simplifications render certain aspects of the problem more amenable to others because they 
are rather insensitive to the assumptions. The scaling of tcoUapse is shown to depend on the behavior 
of a few modes setting the characteristic time scale of the problem - therefore the prediction is not 
affected by the underestimation of the coupling to other modes. 

Even within these simplications, we have obtained good quantitative predictions for the onset of 
instability, the scaling of the waiting time tcoUapse (when the condensate implosion really begins after 
the inversion of the scattering length) with the scattering length, and also for the number of particles 
in a jet as a function of Tevoive, when the interaction between atoms is switched off. 

Another success of the model is to provide a simple explanation for the widely different appearance 
of bursts and jets. As remarked earlier, jets may only appear if the turn - off time T^voive is earlier than 
the formation of the remnant, because once the condensate is stable again, there are no more frozen 
modes to thaw, but, on the other hand, jets will appear for Tevoive < tcoiiapse-, when the condensate 
has not yet shed any particles. Also jets must be less energetic than bursts, since they are composed 
of lower modes. 

Considering the success these simple ideas and light calculations brought about we believe our 
approach might have captured the essence of the physics behind these phenomena. The physical 
paradigms we used in bringing forth these ideas also suggest that understanding the basic mechanism 
of important processes in cosmology, critical dynamics, and Bose-Einstein condensation may share 
more than a superficial ground. 
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6 Appendix A: the classical GPE 



In this appendix we shall consider the purely classical GPE, meaning that we shall disregard any 
back reaction from fluctuations. The conclusion of this analysis is that for short times of the order of 
1ms, it is possible to consider the modulus of the condensate as time independent, and its phase as 
space independent. This observation will lead to a substantial simplification of the equations for the 
fluctuations in this regime. 

Let us begin from the equations 



^ + vGv$o + 1 (v^e) #0 = 

It is suggestive to rewrite these equations in the following way. The operator 



(102) 
(103) 



A = ^ + (ve)v 



(104) 



looks a lot as a material derivative with respect to a fluid flowing at each point with velocity v = V0. 
It seems natural to change from the eulerian coordinates x to Lagrangian coordinates q. That is, for 
given q we define the function x = x {q, t) as the solution to the system 



dx {q, t) 



dt 



= V^e{x{q,t),t), x{q,0)=q 



Observe that 



A 



It is also convenient to define the fields 



r 



d_ 
di 



dx^ 
dq^ 



(105) 



(106) 



(107) 



and their inverses J°. At any given time, the q"" are a set of curvilinear coordinates, with metric 



ds =gapdq°'dq' 



(108) 



The wave function transforms as a scalar under this coordinate change. We also define 

J = det 4 



(109) 
(110) 



Observe that 



22 



dJ_ 

dt 



J J"" 



dq'^dx^ 



The volume element is g = det ga/s = J^, and the Laplacian 



V5 



This suggests writing 



To get the equation 



$(1 



5*0 

dt 



= 



(111) 



(112) 



(113) 



(114) 



Writing the dynamics this way, the role of the phases is hidden in the time dependent curvilinear 
coordinate system, since we get 

*o = *o (?) = ^e-^V2 (115) 

where we are using the fact that initially the condensate corresponds to a noninteracting gas. 
To be definite, consider a situation with Nq = 16000 and a = 30 oq. We obtain 

Nqu = An— No = 104 (116) 
aho 

Since we expect that, at least initially, H^q ~ 3$o/2, the dynamics of the phase near the origin is 
dominated by the u^q term. Concretely, observe that 



de 

'dt 



'dt 



(vey 



So the equation for the phase reads 



Approximate 



Observe that 



Q{q,t) = no-t 



(117) 



(118) 



(119) 



Ve = 2tuq^l (q) 



(120) 



So 
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(121) 



This term is therefore negHgible against u^q ii q ^ 1 or q < {t^/Nou) ^ ; it is small everywhere if 
t < [\/Nou) ^ . This is enough to justify the formal procedures in the main body of this paper. 
The equation which defines the coordinate change from a; to g becomes 



dx {q, t) 



dt 



= 2tuq^l (q) , X {q, 0) = q, 



(122) 



with solution 



X 



1 + fu^l (q) 



So again, we see that during the first millisecond, we may approximate x 
establishes the conclusions anticipated at the beginning of this section. 



(123) 

q everywhere. This 



7 Appendix B: exact solutions of the mode equation in the dynam- 
ical case 

In this Appendix we shall derive closed form solutions for the evolution equations for quantum fluc- 
tuations after tcollapse 



Therefore 



Call 



Therefore 



^7. 



^« - ( ^ ) exp (-i/r) 



dt 



-\- 

dt^ 



- ( ^ ) exp (-t/r) 



Cn = 



dt 2r dC 



(124) 
(125) 

(126) 
(127) 
(128) 



and 



^d 



Eft, — 



auJz\ f C 



a J VCo 



Cn = 



(129) 



so 
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Choose 

Co^ = l^^f^^ (131) 

to find (recall that the c^j are self-adjoint) 

Cn = oinhiTEn (C) + Q;i-f-2jr£;3 (C) (132) 

We see the two basic behaviors discussed in the previous section. For t « r, ( » 1, the Bessel 
functions behave as real exponentials. When i — ^ oo, ^ the modes oscillate with frequency 
±Efi. Modes with Efi > auz/a are never frozen. Lower modes start frozen at t = 0, but thaw when 
C ~ 2tE^, or 

e(-t/r) ^ ^ (133) 



The hfi coefiicients are given by 



1 rjc^ —C dc^ 



The Cfi and 6^ coefficients must be continuous, so 



anhiTEfi (Co) + ail-2iTEa (Co) = (135) 

-Co 



OlH-^hirEn (Co) + Oi\-^I-2iTEn (Co) 

where Cfj and 6ij are the result of the evolution up to t collapse-, namely 



hn (136) 



2i 

Using the Wronskian 



If t \ 1 («n - at) (0) 

CrJ = - (^a^ + al j (0) cosh afitcollapse + JZ ^^^^ (^ntcollapse (137) 

= j: {aft- aVj (0) cosh a^jicoiiapse + ^'&n (ori + at) (0) Sinh. aritcollapse (138) 



/> N r /AN 7- \ d T /J- \ „.sinh27rr£^fj 

^2ir£;3 (Co) -J^I-2iTEn (Co) - -^-2iTE. (Co) -^^^irE^ (Co) = -2? ^j^^ (139) 



we get 



Call 



ivrCo 



2 sinh 2'KTEfi 



-T:I-2iTEa (Co) Cf; + I-2iTEfi (Co) — r^^'n 

"C Co 



(140) 
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Then 



^^(C,Co) = 
w^i(C,Co) = 



2 sinh 2'ktE, 



hirEfi (C) I-2iTEfi (Co) - I-2iTE^ (C) hirEfi (Co) 



So, writing 



Then 



2 sinh 2'KTEfj 



hirEfi (C) -^I-2irEfi (Co) - I-2irEfi (C) J^hirEfi (Co) 

w^2(C,Co) = ^w^(C,Co) 
W3(C,Co) = ^w-i(C,Co) 



Cn = H^i (C, Co) C^J + — (C, Co) bn 
Co 



-c 

2E^jr 



2E^T 

W^3(C,Co)cs + — ^VF2(C,Co)6fi 

Co 



Cn = ^ (an + 4) ; = ^ ("^ ~ 4) 



2E^T 

an it) = Wi iC, Co) CH + —^WiC, Co) bn 

Co 



iC 



2EnT 



2E^T - ' 

(C, Co) Cn + —^VF2(C, Co) 

Co 



an = fnan (0) + gnai (0) 



an(i) = /n(i)an(0)+5n(i)4(0) 



hit) = 2^i(C>Co)(/n + ffi) 

-^W^(C,Co) (/n-5i) 

Co 

--^1^3 (C, Co) ifn + rn) 

-l^W2iC,Co) ifn-9*n) (152) 
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9n{t) = lWl{C,Co){fH+9n) 

VF3(C,Co)(/^+5n) 



AEht 

+^VF2(C,Co)(/5-5n) (153) 



This result is used to build the plot in Fig. 6. 
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8 Figure Captions 



8.0.1 Fig. 1 

Effective scattering lengtfi as a function of applied magnetic field (unlike in the body of the paper, we 
follow here the usual sign convention) . The scattering lenght is measured in multiples of Bohr radius 
oq = 0.529 lO^^^'m. The magnetic field is measured in Gauss. In the Bose Novae experiment 0121) 
the condensate is prepared at Bzero and then evolved to a larger field. In later experiments [51 [S] the 
field was turned to a lower value, close to Bp^ak for a short time, and somewhere between Bp^ak and 
Bzero for tlie duration of the experiment. 

8.0.2 Fig. 2 

Qualitative evolution of the effective scattering length as a function of time in the Bose Novae experi- 
ment 012] (unlike in the body of the paper, we follow here the usual sign convention). The scattering 
lenght in measured in multiples of Bohr radius bq = 0.529 10~^''m, time is measured in ms. The 
condensate is prepared at a = 0, and then the effective interaction is made attractive. After a time 
Tevoivei the interaction is made repulsive, which allows the condensate to expand and aids visualization. 

8.0.3 Fig. 3 

Plot of tcoUapse (in nis) against acoiiapse (in multiples of Bohr radius oq = 0.529 10~^°m). We plot the 
scaling law Eq. (|56|) (full line) and compare it against the experimental data for Nq = 6000 as reported 
in Refs. [2| (small black dots), the tjyL ~ {uNo)~^ prediction (suitably scaled) as given in 
(dashed line) and the results of numerical simulations reported in |2J| (large grey dots). While all 
three theoretical predictions may be considered satisfactory, the tNL ~ {uNo)~^ fails to describe the 
divergence of tcoUapse as the critical point is approached. The results of numerical simulations based 
on an improved Gross-Pitaievskii equation tend to be systematically above the experimental results. 
In a classical instability, the unstable modes must grow from zero, while in a quantum instability, 
they are always seeded by their own zero - point fluctuations, which speeds up the development of 
the instability. Therefore, the fact that numerical simulations tend to overestimate tcoUapse may be a 
further indication of the quantum origin of the phenomenon. 

8.0.4 Fig. 4 

The evolution for the number of particles in the condensate N (t) as a function of time (measured 
in ms) assumed for the calculation of the number of particles in a jet. For comparison, we have 
superimposed the data from Fig. lb of fll. The agreement is satisfactory for our purposes, as we shall 
not consider the evolution beyond Tevoive = 12ms nor the formation of a remnant. 

8.0.5 Fig 5 

A schematic depiction of the evolution of the proper frequencies of different modes in time, both in 
arbitrary units. Short wavelength modes remain stable throughout, but the sudden change in proper 
frequency when the interaction is switched on induces particle creation in them. These particles are 
perceived as bursts. Long wavelength modes actually become unstable (or frozen) until the condensate 
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has lost enough atoms. During this period they are amphfied. When they become stable again, they 
are seen as a secondary emission, or jet. 

8.0.6 Fig. 6 

The evolution of the number of particles in a jet as a function of the time t^voive (measured in ms), 
for No = 16,000, uJradiai = 110 Hz, ujaxiai = 42.7 Hz, a = 36oo, and k = 0.46 (see Eq. ^). The 
evolution of N {t) in time is depicted in Fig. 4. To provide a visual reference we have superimposed 
the data in fig. 6 of 1 . The discrepancy between theory and experiment beyond tevoive = 6ms may 
be attributed to an overestimation of the condensate-noncondensate coupling in neglecting the change 
in the shape of the condensate. 
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